Daphnias: from the individual based model to the large population equation

The class of deterministic 'Daphnia' models treated by Diekmann et al. (J Math Biol 61: 277-318, 2010) has a long history going back to Nisbet and Gurney (Theor Pop Biol 23: 114-135, 1983) and Diekmann et al. (Nieuw Archief voor Wiskunde 4: 82-109, 1984). In this note, we formulate the individual based models (IBM) supposedly underlying those deterministic models. The models treat the interaction between a general size-structured consumer population ('Daphnia') and an unstructured resource ('algae'). The discrete, size and age-structured Daphnia population changes through births and deaths of its individuals and throught their aging and growth. The birth and death rates depend on the sizes of the individuals and on the concentration of the algae. The latter is supposed to be a continuous variable with a deterministic dynamics that depends on the Daphnia population. In this model setting we prove that when the Daphnia population is large, the stochastic differential equation describing the IBM can be approximated by the delay equation featured in (Diekmann et al., l.c.)

Numerical equilibrium analysis for structured consumer resource models

In this paper, we present methods for a numerical equilibrium and stability analysis for models of a size structured population competing for an unstructured resource. We concentrate on cases where two model parameters are free, and thus existence boundaries for equilibria and stability boundaries can be defined in the (two-parameter) plane. We numerically trace these implicitly defined curves using alternatingly tangent prediction and Newton correction. Evaluation of the maps defining the curves involves integration over individual size and individual survival probability (and their derivatives) as functions of individual age. Such ingredients are often defined as solutions of ODE, i.e., in general only implicitly. In our case, the right-hand sides of these ODE feature discontinuities that are caused by an abrupt change of behavior at the size where juveniles are assumed to turn adult. So, we combine the numerical solution of these ODE with curve tracing methods. We have implemented the algorithms for “Daphnia consuming algae” models in C-code. The results obtained by way of this implementation are shown in the form of graphs

Connectivity, neutral theories and the assessment of species vulnerability to global change in temperate estuaries

One of the main adaptation strategies to global change scenarios, aiming to preserve ecosystem functioning and biodiversity, is to maximise ecosystem resilience. The resilience of a species metapopulation can be improved by facilitating connectivity between local populations, which will prevent demographic stochasticity and inbreeding. The objective of this investigation is to estimate the degree of connectivity among estuarine species along the north-eastern Iberian coast, in order to assess community vulnerability to global change scenarios. To address this objective, two connectivity proxy types have been used based upon genetic and ecological drift processes: 1) DNA markers for the bivalve cockle (Cerastoderma edule) and seagrass Zostera noltei, and 2) the decrease in the number of species shared between two sites with geographic distance; neutral biodiversity theory predicts that dispersal limitation modulates this decrease, and this has been explored in estuarine plants and macroinvertebrates. Results indicate dispersal limitation for both saltmarsh plants and seagrass beds community and Z. noltei populations; this suggests they are especially vulnerable to expected climate changes on their habitats. In contrast, unstructured spatial pattern found in macroinvertebrate communities and in C. edule genetic populations in the area suggests that estuarine soft-bottom macroinvertebrates with planktonic larval dispersal strategies may have a high resilience capacity to moderate changes within their habitats. Our findings can help environmental managers to prioritise the most vulnerable species and habitats to be restored

Effects of aging and links removal on epidemic dynamics in scale-free networks

We study the combined effects of aging and links removal on epidemic dynamics in the Barab\'{a}si-Albert scale-free networks. The epidemic is described by a susceptible-infected-refractory (SIR) model. The aging effect of a node introduced at time $t_{i}$ is described by an aging factor of the form $(t-t_{i})^{-\beta}$ in the probability of being connected to newly added nodes in a growing network under the preferential attachment scheme based on popularity of the existing nodes. SIR dynamics is studied in networks with a fraction $1-p$ of the links removed. Extensive numerical simulations reveal that there exists a threshold $p_{c}$ such that for $p \geq p_{c}$, epidemic breaks out in the network. For $p < p_{c}$, only a local spread results. The dependence of $p_{c}$ on $\beta$ is studied in detail. The function $p_{c}(\beta)$ separates the space formed by $\beta$ and $p$ into regions corresponding to local and global spreads, respectively.Comment: 8 pages, 3 figures, revtex, corrected Ref.[11

On models of physiologically structured populations and their reduction to ordinary differential equations

Considering the environmental condition as a given function of time, we formulate a physiologically structured population model as a linear non-autonomous integral equation for the, in general distributed, population level birth rate. We take this renewal equation as the starting point for addressing the following question: When does a physiologically structured population model allow reduction to an ODE without loss of relevant information? We formulate a precise condition for models in which the state of individuals changes deterministically, that is, according to an ODE. Specialising to a one-dimensional individual state, like size, we present various sufficient conditions in terms of individual growth-, death-, and reproduction rates, giving special attention to cell fission into two equal parts and to the catalogue derived in an other paper of ours (submitted). We also show how to derive an ODE system describing the asymptotic large time behaviour of the population when growth, death and reproduction all depend on the environmental condition through a common factor (so for a very strict form of physiological age).Peer reviewe

Finite dimensional state representation of physiologically structured populations

In a physiologically structured population model (PSPM) individuals are characterised by continuous variables, like age and size, collectively called their i-state. The world in which these individuals live is characterised by another set of variables, collectively called the environmental condition. The model consists of submodels for (i) the dynamics of the i-state, e.g. growth and maturation, (ii) survival, (iii) reproduction, with the relevant rates described as a function of (i-state, environmental condition), (iv) functions of (i-state, environmental condition), like biomass or feeding rate, that integrated over the i-state distribution together produce the output of the population model. When the environmental condition is treated as a given function of time (input), the population model becomes linear in the state. Density dependence and interaction with other populations is captured by feedback via a shared environment, i.e., by letting the environmental condition be influenced by the populations' outputs. This yields a systematic methodology for formulating community models by coupling nonlinear input-output relations defined by state-linear population models. For some combinations of submodels an (infinite dimensional) PSPM can without loss of relevant information be replaced by a finite dimensional ODE. We then call the model ODE-reducible. The present paper provides (a) a test for checking whether a PSPM is ODE reducible, and (b) a catalogue of all possible ODE-reducible models given certain restrictions, to wit: (i) the i-state dynamics is deterministic, (ii) the i-state space is one-dimensional, (iii) the birth rate can be written as a finite sum of environment-dependent distributions over the birth states weighted by environment independent 'population outputs'. So under these restrictions our conditions for ODE-reducibility are not only sufficient but in fact necessary. Restriction (iii) has the desirable effect that it guarantees that the population trajectories are after a while fully determined by the solution of the ODE so that the latter gives a complete picture of the dynamics of the population and not just of its outputs.Peer reviewe